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The fluid flow and forced convective heat transfer characteristics of 
a uniformly heated right circular cylinder resting on a flat plate in a 
crossflow of air were experimentally investigated. 

The circumferential variation of the Nusselt number on the cylinder 
surface was obtained at Reynolds numbers ranging from 53,000 to 153,000 
using a liquid crystal thermographic technique. Pressure coefficients 
were obtained at five-degree intervals around the cylinder circumference 
at Reynolds numbers of 105,000 and 156,000. Flow visualization was con- 
ducted at a Reynolds momber of 10,000 in order to qualitatively examine 
the flow field. 

A qualitative theory based upon a compilation of previous results 
reported for flow around a "free" cylinder and flow over forward and 
rear-facing steps was developed. 

Experimental results obtained were in agreement with the proposed 
theory, and heat transfer and pressure data correlated well with one 
another. In the region between 0 and 180 degrees, the heat transfer and 
pressure results correlated quantitatively with the results obtained by 
Meyer for a "free" cylinder at Reynolds numbers above 300,000. 
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I . INTRODUCriON 



The primary purpose of this work was to investigate the forced convec- 
tive heat transfer and fluid flow characteristics of a uniformly heated 
right circular cylinder resting on a flat plate in a crossflow of air. 

This particular flow situation models several situations of practical 
interest/ notably, the effect of wind on rocket motor storage containers 
lying on the ground. In this regard, this investigation supports the 
previous work reported by Wirzburger [1]*. 

Local Nusselt numbers were obtained at various circumferential loca- 
tions on the cylinder for Reynolds numbers in the range of 53,000 to 153,000 
using the liquid crystal technique developed by Field [2] . A four-inch 
diameter cylinder was rigidly attached to a flat plate and placed in a lew 
speed wind tunnel with a turbulence intensity of approximately 0.5 to 0.7 
percent. The test cylinder was constructed of a carbon-impregnated resis- 
tive paper having an electrical resistivity of approximately one ohm-inch. 
The cylinder was coated with microencapsulated cholesteric liquid crystals 
and packed inside with glass wool for the purposes of both minimizing 
conduction and internal free convection losses, and providing resistance 
to deformation. The cylinder surface was uniformly heated by passing an 
electrical current through the resistive paper, and data v;ere obtained 
by observing the position of the isotherms indicated by the liquid crystals. 

Pressure coefficients were also obtained on the surface of the cylinder. 
A rotatable four- inch diameter aluminum cylinder was fitted with four 0.025 
inch static pressure taps 90 degrees apart. The cylinder was sealed to 
the flat plate by moans of a Teflon which allowed the c:yixndor to 

♦Numbers in brackets refer to references at end of thesis. 
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rotate, but did not allow flow to pass between the plate and the cylinder. 
Pressure data were obtained using a referenced against 
atmospheric pressure. Stagnation pressure and test section static pressure 
were obtained from existing pressure taps in the wind tunnel which were 
connected to the same manometer bank. 

Flow visualization was conducted in a water tunnel using a 1.5-inch 
diameter plexiglas cylinder sealed against an aluminum flat plate. The 
system was scaled for dynamic similitude with the heat transfer experimental 
apparatus. The cylinder was equipped with three dye ports, two of which 
were located approximately 30 degrees either side of the attachment point, 
and the other located at approximately the 90 degree mark. Due to limita- 
tions in the water pump, the maximum Reynolds number which could be 
attained with this apparatus was approximately 10,000. 

A qualitative theory of the flow field and accompanying heat transfer 
was proposed based upon previous investigations of flow around a "free*' 
cylinder combined with investigations of flow over forward and rear-facing 
steps. 

The heat transfer and pressure data obtained were in excellent agree- 
ment with the proposed theory, and the pressure data correlated well with 
the heat transfer data. The flow visualization showed agreement with 
portions of the proposed theory, but did not support the heat transfer 
and pressure data on certain regions of the cylinder. This is most 
probably because of the wide difference in the Reynolds number between 
these experiments. 
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11 . THEORY 

A thorough search of the existing literature indicated that this 
particular flow situation had not been examined. The study of flow around 
a ”free” cylinder has been studied extensively^ and an excellent bibliog- 
raphy on the subject has been presented by Meyer [3] . Subsequent to 
Meyer's work. Field [2] investigated both subcritical and critical flow 
around a right circular cylinder in crossflow using the liquid crystal 
technique employed in this investigation. A brief description of these 
two flow situations is presented below. 

In subcritical flow, a laminar boundary layer develops at the forward 
stagnation point of the cylinder and grows to a maximum thickness at an 
angular location of 80 to 85 degrees. Concurrently, the heat transfer 
coefficient decreases as the boundary layer thickness increases in this 
region. At 80 to 85 degrees, the kinetic energy of the moving fluid is 
no longer sufficient to overcome the adverse pressure gradient present 
on the cylinder's surface, and the laminar boundary layer separates. 

Upon separation, a wake is formed in the rear portion of the cylinder, and 
the turbulent action in the wake tends to circulate cooler fluid to the 
surface of the cylinder. As a result of this action, the heat transfer 
coefficient increases in this area. 

In critical flow, which occurs at Reynolds numbers above approximately 
120,000, the laminar boundary layer again grows from the forward stagnation 
point to the 80 to 85 degree point where separation again occurs. However, 
in this case the kinetic energy of the fluid is enough to overcome the 
adv^erse pressure gradient, and the flow reattaches itself to the cylinder 
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surface at a point approximately 15 degrees aft of the laminar separation 
point. This separation highly turbulent 

mixing and results in a sudden increase in heat transfer. From the point 
of reattachment, a turbulent boundary layer begins to grow to a point in 
the region of 110 to 130 degrees where it separates. The heat transfer 
coefficient decreases as the turbulent boundary layer thickens. Upon 
separation of the turbulent boundary layer, the action of the wake causes 
an increase in the heat transfer coefficient. 

It should be noted that the above description is confined to a right 
circular cylinder which is freely suspended in crossflow such that flow 
can pass over both sides. Hence, the flow pattern is spacially symmetrical. 
There was no reason to expect this to be the case in this investigation, 
however, and an attempt was made to gain some understanding of what results 
might be expected in the regions of the cylinder near the flat plate. This 
was accomplished by examining the cases of flow over forward-facing and 
rear-facing steps. 

The problem of flow over forward-facing and rear-facing steps has been 
examined in various lights. Several investigators have focused on the 
determination of the location of the separation points and reattachment 
points [4-9] • Others have examined the heat transfer and pressure dis- 
tribution within the separated regions [10-17] . However, these investiga- 
tions were primarily concerned with the regions of the flat plate either 
forward or aft of the step and did not deal with the step itself. Since a 
cylinder resting on flat surface represents a type of "step", this region 
was of most interest in this work. A few investigators did conduct flow 
visualization experiments in conjunction with their other goals, however, 
and their results proved to be of benefit in predicting the nature of the 
flow in the regions of the cylinder near the flat plate. 
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Of particular interest were the findings of Abbott and Kline [5] for 
rear-facing steps, and those of [10] for forward- 

facing steps. The latter investigated the zone of turbulent boundary layer 
separation ahead of a forward-facing step at Reynolds numbers of 17,000 to 
420,000 based upon step height. Although their visualization technique was 
not described, their experiments were conducted in air. The results of 
this investigation are shown in Figure 1. It should be noted that, in 
addition to the normal separation vortex (points 1-2-1) , another small 
"trapped" vortex was noted in the region between points 3 and 4. 

Abbott and Kline conducted flow visualization experiments using the dye 
injection method on a variable height backward-facing step placed in a 
water tunnel. The boundary layer at separation from the step was turbulent 
as is the case for critical flow around a cylinder. The experiments were 
conducted in a Reynolds number range of 20,000 to 50,000 based upon the 
height of the water channel above the step. Their results indicated that 
three regions of flow existed downstream of the step. Region I was 
characterized by one or more vortices rotating about an axis normal to the 
vertical wall and parallel to the floor. These vortices were three- 
dimensional in nature, and adjacent vortices were counter-rotating and not 
necessarily the same size. Region II contained the classical pattern of 
separation, i.e. flow moving two-dimensionally upstream next to the wall 
and dov;nstream adjacent to the through flow. Region III was of no interest 
to the present investigation because it did not directly affect the step. 
The three regions are depicted both in plan and elevation views in Figure 
2 . 

From a compilation of the results described above and those observed 
for a "free" cylinder, the following theory was developed. 
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Figure 1. Flow pattern ahead of a step as observed by 
Luzhanskiy and Solntsev. 
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Figure 2. Flow pattern at the rear of a step as observed 
by Abbot: t and Kline. 
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As the flow proceeds along the flat plate, a boundary layer will develop 
and grow on the plate until, at some point upstream of the cylinder, this 
boundary layer will separate. At the flow velocities to be used in the 
experimental portion of this investigation, this boundary layer can be 
anticipated to be laminar. From this point, it is theorized that the 
flow will be characterized by a behavior similar to that observed by 
Luzhanskiy and Solntsev, i.e. a portion of the flow will reattach to the 
cylinder at some point on the forward portion of its circumference, and 
the remainder of the flow will assume the classical separation pattern with 
the addition of a system of "trapped" vorticies near the intersection of 
the cylinder and the flat plate. The flow which reattaches itself to the 
cylinder then will assume the classical pattern described previously for 
flow around a "free" cylinder, with the differentiation between subcritical 
and critical flow being dependent upon the Reynolds number. Regardless of 
whether the flow is subcritical or critical, it will separate. Upon 
separation, it is assumed that a pattern similar to that observed by 
Abbott and Kline will be present. 

From the above theory of the flow phenomenon, the heat transfer results 
can be predicted. For ease of explanation the cylinder will be described 
as consisting of three regions. Region I consists of that portion of the 
cylinder circumference betv/een the point of attachment of the cylinder to 
the flat plate and the point at which the separated flow from the plate 
reattaches to the cylinder. This latter point will hereafter be referred 
to as the stagnation point for simplicity. Region II consists of that 
portion of the cylinder circumference between the stagnation point and the 
point of final boundary layer separation. Region III consists of that 
portion of the cylinder circumference aft of the final sei)aration point. 
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The phenomenon which is proposed to occur in Region I is depicted in 
Figure 3. Since all the fluid zone is "trapped," 

heat liberated from the cylinder surface must be conducted from vortex to 
vortex within the separation region and ultimately conducted to the fr6e 
stream through the free shear layer separating the separation zone from 
the free stream. This being the case, the bulk temperature in the larger 
vortices must be less than the bulk temperature in the smaller vortices 
in order for a driving potential for conduction to exist. Convection, 
however, also plays a role in this process in that a redistribution of 
energy within individual vortices is accomplished by convection. It is 
assumed that the vortices will become smaller and weaker as the attachment 
point is approached such that, very near this point, the heat transfer will 
approach a pure conduction case. Hence, the resistance to heat transfer 
should increase from the stagnation point to the attachment point. If the 
thermal resistance increases, the wall temperature must also increase in 
order to insure that a constant heat flux is transferred to the fluid. If 
the wall temperature increases and the free stream temperature and surface 
heat flux remain constant, the heat transfer coefficient must decrease. 
Hence, in the absence of other effects, the heat transfer coefficient 
should decrease from a relative maximum at the stagnation point to a zero 
value at the attachment point (assuming the flat plate is adiabatic) . 
Superimposed upon on this phenomenon, hov;ever, is the behavior of the heat 
transfer coefficient near the intersection of two vortices. The region 
between points 1 and 3 of Figure 3 is illustrative of this behavior. At 
point 1, the stagnation point, the fluid in the vortex which impacts on 
the cylinder surface is relatively cool because some of its energy has 
just been conducted to the free stream. As the fluid proceeds toward 



15 



m 




16 



point 2 , it is heated and its thermal resistance is increased; hence, the 
heat transfer coefficient is decreased' Since adjacent 
vortices are counter-rotating, a similar behavior should occur between 
points 3 and 2. Therefore, point 2 should appear as a relative minimum 
in the heat transfer coefficient, and the heat transfer coefficient should 
then increase to a relative maximum at point 3. At point 3, the cycle 
described above should begin again, this time involving points 3, 4, and 5. 
It should be noted that, because the vortices become smaller and weaker 
as the point of attachment is approached, this secondary effect may become 
experimentally imperceptible in the region near this point. It should also 
be noted that the number of vortices depicted in the schematic does not 
consititute a prediction of how many will actually be present. 

Region II is characterized by the results reported by Field [2] for 
subcritical or critical flow, whichever may be the case. 

Region III is theorized to indicate an initial increase in the heat 
transfer coefficient due to the turbulent action of the wake. However, 
at some point, this increase will be overridden by the phenomenon described 
for Region I, and the heat transfer coefficient will decrease to a zero 
value at the attachment point. 

Schematics depicting the theorized flow situation and accompanying heat 
transfer data for both subcritical and critical flows are shown as Figures 
4 and 5 respectively. Tlie locations at which various phenomenon occur 
are conjectural and should not be construed as being exact predicted 
values . 
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Ficjurr^ A, Proposed flow pattern for :;ubcr i Li cal Reynolds number 
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Ill . EXPERIMENTAL APPARATUS 



A. HEAT TRANSFER EXPERIMENT 

1. Liquid Crystals 

Microencapsulated cholesteric liquid crystals were used in this 

investigation. This substance/ when applied to a dark background/ provides 

a visual display of temperature by means of a change of color. The surface 

coated with liquid crystals traverses the visible color spectrum within a 

o 

certain temperature band with an accuracy of approximately 0.2 F. The 

size of this band is determined by the formulation of the cholesterol esters 

in the mixture. Only one liquid crystal formulation was used in this 

investigation/ NCR designation S-43. The batch used had been calibrated 

previously by Field [2], and the results of this calibration indicated 

o 

that liquid crystal S-43 transitioned to a red color at 109.8 F/ to a 

o o 

green color at 111.0 F, and to a blue color at 112.3 F. An excellent list 
of references on the theory and use of liquid crystals can be found in 
Field's work. 

2 . Cylinder Design and Construction 

The cylinder utilized in the collection of heat transfer data was 
similar in design to that used by Field. The actual working section was 
constructed from a 12-inch square sheet of 7\rmstrong Temshect, a carbon- 
impregnated, electrically resistive paper having an electrical resistance 
of approximately 25 ohms per square and a nominal thickness of 0.039 inches. 
The heat generated when a constant current is passed through the paper is 
uniform to within two percent from point to point over a large area. This 
sheet was fashioned into a cylinder having an outside diameter of 4.03 
inches and a height of 12 inches. 
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Alumimam conducting tape was applied as electrodes circumferentially 
on the inner surface of this cylinder to insure uniform application of volt- 
age, These electrodes were located three inches from the top and bottom of 
the cylinder in order that edge conduction effects would be eliminated. 

This established a six inch height of cylinder in which data would be taken. 
Silver conducting paint was applied to the interface between the tape and 
the paper to insure intimate electrical contact was maintained between the 
two surfaces. 

The base assembly consisted of two 4.03-inch diameter wooden 
cylinders, each 14.75 inches long, rigidly attached to each other by an 
eight-inch long wooden bar. The outside surface of this bar was of the 
same curvature as the cylinders in order to insure that no discontinuities 
existed near the point of attachment to the flat plate. The purpose of 
this bar was to maintain alignment between the two base cylinders, and to 
allow the Temsheet cylinder to be taped along its length to the bar. The 
entire assembly was slotted to accept a 1/4-inch square plexiglas rod along 
its length. This rod accepted screws from the flat plate to insure that 
an acceptable seal was maintained at the point of attachm^ent. The base 
assembly is shov,Ti in Figure 6. 

Each base cylinder v/as of sufficient length such that it would 
extend through the floor/ceiling of the v;ind tunnel where it could be 
clamped thus providing added rigidity to the assembly. The other ends of 
the base cylinders were turned down approximately 0.04 inches over a two- 
inch length. This allowed the paper to be attached to the bases with 
double-backed tape and present a constant diameter of 4.03 inches along 
the entire length of the assembly. The bottom base was drilled through 
its axis v;ith a one^inch diameter hole through which the electrica.1 leads 



21 













Figure 6. Temsheet cylinder base assembly 
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could be passed to the power supply; and the top base was drilled through 
its axis with a two-inch diameter hole to provide a means of filling the 
Temsheet cylinder with glass wool. The glass wool served the dual purpose 
of providing insulation against possible internal heat losses due to con- 
duction or natural convection and of providing the paper cylinder with 
added resistance to deformation. 

Two coats of slurry-based liquid crystal S-43 were applied to the 
Temsheet cylinder, and it was marked at five-degree intervals for accurate 
location of the liquid crystal isotherms during the experiment. 

The flat plate was constructed from a sheet of 1/2-inch plexiglas 
and was 36 inches long and 32 inches high. Plexiglas was selected for this 
constructipn in order that the entire circumference of the cylinder could 
be viewed from outside the wind tunnel. The leading edge of the plate was 
beveled approximately 14 degrees to insure that a sharp leading edge v;as 
presented to the flow. The point of attachment of the cylinder was located 
on a vertical line 18 inches downstream of the leading edge. At this 
point, the plate was notched along its height to accept the plexiglas red 
on the cylinder assembly, and attachment was accomplished by four screws. 
This design insured that flow could not pass between the cylinder and the 
plate. 

The wind tunnel in which this cylinder assembly was placed was the 
low speed Aerolab wind tunnel located in the Aeronautics Laboratory in 
Halligan Hall, U. S. Naval Postgraduate School. This tunnel is powered by 
a 100 hp electric motor and is capable of producing wind velocities of up 
to 200 mph v;ith a clear test section. The tunnel has a turbulence level 
as measured by hot wire anemometer of approximately 0.5 to 0.7 percent. A 
schematic of the wind tunnel v;ith pertinent dimensions is shown in Figure 7. 
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Fi jure 7. Scheinatic diagram of the wind tunnel. 




The combined cylinder and plate assembly was located in the wind 
tunnel test section such that the longitudinal axis of the flat plate lay 
along a line eight inches from one wall of the test section. By off- 
^ centering the assembly in this manner, the effects of the flow along the 
wind tunnel wall on the cylinder were effectively eliminated. 

B . PRESSURE EXPERIMENT 

In order to. obtain static pressure measurements of this particular 
flow situation, a 4.03-inch diameter cylinder was constructed of aluminum. 
This cylinder was fitted with four 0.025-inch surface pressure taps spaced 
at 90-degree intervals around the circumference of the cylinder at approxi- 
mately mid-height of the wind tunnel test section. In order that this 
cylinder could be rotated and a seal still be maintained, a Teflon strip 
was inserted into the notch in the flat plate described previously. This 
strip provided for a snug fit between the cylinder and the plate, and the 
flexibility of the Teflon allowed for relatively easy rotation. The 
cylinder was scribed as five-degree intervals about its circumference at 
the point where it intersected the wind tunnel test section floor. 

Reference lines were scribed on the floor at 90-degree intervals so that 
the matching of the scribe marks on the cylinder with these lines provided 
a means of accurately establishing the angular location of the pressure 
taps. The pressure cylinder is shov/n in Figure 8. 

The surface pressure taps were connected to a water manometer bank by 
means of plastic tubing. The monometer bank was referenced against 
atmospheric pressure. Wind tunnel plenum chamber static pressure and test 
section static pressure were measured on the same manometer bank through 
pre-existing connections on the wind tunnel. 
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G. FLOW VISUALIZATION EXPERIMENT 



The apparatus used for flow visualization consisted of a converging 
wooden water channel 15 feet in length and 13 inches deep with a plexiglas 
test section 13 inches long, 12 inches wide, and 13 inches deep, A hollow, 
1.5-inch diameter plexiglas cylinder was attached between the test section 
walls with its longitudinal axis located eight inches from the test section 
floor. In order to preserve dynamic similitude, a 1/4-inch thick steel 
flat plate was constructed to the scale of the ratio of the water tunnel 
cylinder diameter to the wind tunnel cylinder diameter (0,375:1), Thus, 
this plate was 13,5 inches long and 12 inches wide. The leading edge was 
beveled to a 14 degree angle as was the larger flat plate used in the wind 
txinnel. The flat plate was placed atop the cylinder, and the point of 
attachment was sealed by a layer of tape. The cylinder was fitted with 
three dye ports which were connected by rubber tubing to three separate 
dye reservoirs with individual valving located approximately three feet 
higher than the cylinder. This additional height insured that a constant 
head was maintained on the injection system, and the individual valving was 
necessary to allow each port to be operated separately in order that the 
flow field at various sections of the cylinder could be observed independently. 

V7ater was pumped into the channel by means of an irrigation pump directly 
coupled to a VJestinghouse Type SK dynamometer running as a motor. V^ith the 
motor running at 1500 RPM, the pump output was approximately 1000 gallons 
per minute. The depth of water in the channel was controlled by a tilting 
weir at the end of the channel. For this particular experiment, the water 
depth was maintained at 10 inches. This allowed the entire apparatus to 
be immersed to a minimum depth in order that free surface effects could be 
minimized . 
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IV. EXPERIMENTAL PROCEDURES 



A. HEAT TRANSFER EXPERIMENT 

Prior to conducting the experiment, the thermocouple used to measure 
airstream temperature and its ice bath reference junction were examined, 
and the micromanometer used to measure wind velocity was leveled and care- 
fully zeroed. The test cylinder power leads were then connected to a 
Weston ohmmeter, and the cylinder resistance was determined. The leads 
were then switched to a Lambda regulated power supply model LK345A, and 
voltage was applied until the temperature of the entire cylinder surface 
fell within the event temperature range of liquid crystal S-43 causing 
the cylinder to undergo a color change. This procedure was necessary to 
insure that the glass v/ool packing v:as adjusted properly so that regions 
of temperature discontinuity did not exist. The wind tunnel was then 
started and brought to a predetermined velocity at which it was maintained 
throughout the run. Applied voltage was adjusted until distinct positions 
of temperature transition could be observed. After steady state was 
reached, the angular locations and magnitudes of these isotherms were 
recorded in addition to the values of applied voltage and airstream tempera- 
ture. A steady state condition was usually attained in from five to ten 
minutes dependent upon the wind velocity. The voltage was then changed 
until the angular location of the isotherms had changed a significant 
amount, and data were recorded again. This process was repeated until 
isotherms had been observed at reasonable intervals around the circumference 
of the cylinder. Through this procedure, a plot of local Nusselt number 
versus angular location could be quickly and easily obtained. 
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Runs were made at nine Reynolds numbers ranging between 53,000 and 
153,000. Each r\in was accomplished in from 60 to 90 minutes, this time 
being dependent upon Reynolds number. For all runs the points of tempera- 
ture transition were easily distinguishable around the entire circumference 
of the cylinder except in the region from 0 to 45 degrees and in the region 
from 180 to 210 degrees. In these regions, the temperature gradient was 
very shallow, and the entire region appeared isothermal. In addition, 
data were not obtainable in the regions 15 degrees either side of the 
attachment point (the 270 degree mark) because a break in the resistive 
paper was necessary in these regions to facilitate the attachment of the 
cylinder to the flat plate. 

B. PRESSURE EXPERIMENT 

Pressure data v/ere obtained at Reynolds numbers of 105,000 and 156,000. 
Initially, the manometer bank was leveled, all leads were checked for 
tightness, and the air temperature was recorded for evaluation of transport 
properties. The pressure taps were initially located at 0, 90, 180, and 
270 degrees. The wind tunnel was started, and the velocity was increased 
to a predetermined value. Values were recorded for each pressure tap on 
the cylinder, the wind tunnel plenum chamber static pressure, the test 
section static pressure, and atmospheric pressure. The cylinder was then 
rotated in five-degree increments, with data being recorded at each 
increment, until the entire circumference of the cylinder had been traversed. 
Data were then reduced, and a preliminary plot of pressure coefficient 
versus angular location v;as made. This plot v;as examined, and regions 
v/arranting closer investigation v/ere identified. Data were then obtained 
at 2 1/2-degrce increments within these regions. 
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C* FLOW VISUALIZATION EXPERIMENT 



The motor was started and brought up to a speed of 1500 RPM. The 
tilting weir was adjusted such that the height of the water in the channel 
was 10 inches, and the flow was allowed to steady out. The dye was then 
injected through the existing ports in a trial and error fashion until the 
flow could be visualized clearly. Most regions of interest could be 
observed in this manner. In those regions which required closer observa- 
tion, dye was injected directly into the flow by means of a hypodermic 
needle . 

It became necessary to interrupt the procedure periodically due to the 
presence of large air bubbles which coalesced on the flat plate downstream 
of the cylinder. These bubbles were thought to have been caused by 
entrained air in the water which left solution in regions of low pressure. 
These bubbles were swept away by hand before the experiment was continued. 
It should be noted that the effect of these air bubbles was not completely 
deleterious. In certain regions on the cylinder surface where the dye did 
not provide adequate information regarding the existing flow conditions, 
the bubbles could be used as a flow visualization aid. The location of 



the laminar separation point was found in this manner. 
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V. RESULTS 



\ 



A. HEAT TRANSFER EXPERIMENT 

Raw data and sample calculations from the heat transfer experiments 
are contained in Appendix A. This section contains graphical results and 
comparisons of the local Nusselt number versus angular location for various 
Reynolds numbers, and of the average Nusselt number versus Reynolds number. 

Figure 9 shows the results of the heat transfer experiment for Reynolds 
numbers of 53,000, 103,500, and 153,000. The ratio of boundary layer 
thickness on the flat plate to cylinder diameter was approximately 0.3. 

Since the shape of these curves did not vary throughout the Reynolds num- 
ber range tested, the curves presented are considered representative. 

There are several points of interest to be noted from these curves. 

First, the flow around the cylinder exhibited critical behavior at all 
Reynolds numbers tested. Otherwise, there was good agreement between 
these results and the proposed theory. The stagnation point was located 
in the region between 5 and 20 degrees and appeared to move toward 0 degrees 
with increasing Reynolds number. The point of laminar separation was 
relatively stationary at approximately 75 degrees. Reattachment occurred 
consistently at approximately 92 degrees, and final separation occurred 
consistently at approximately 140 degrees. As Reynolds number increased, 
the point of maximum heat transfer shifted from the stagnation point to 
the point of reattachment of the turbulent boundary layer. Of particular 
note are the tv;o points of relative maxima located at approximately 240 
and 300 degrees respectively. Since the point of attachment of the 
cylinder to the flat plate was taken to bo the 270 degree mark, tiiesc 
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Figure 9. iransfer results at Reynolds numbers of 53,000, 103,500, and 153,000. 
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points lie in the regions in which this type of behavior was postulated 
by theory. 

Figures 10 and 11 are a comparison of the heat transfer data obtained 
in this investigation with that obtained by Field [2] at two nearly equal 
Reynolds numbers, one in the subcritical range and one in the critical 
range. Figure 11 also contains a plot of the data obtained by Meyer [3] 
at a Reynolds number of 495,000. Additionally, a comparison of average 
Nusselt number versus Reynolds number for this investigation and Field's 
investigation is shown as Figure 12, The average Nusselt number was 
computed from the formula 




which was numerically integrated by computer using Simpson's Rule of 
Quadrature with a step size of 5 degrees. 

It is readily apparent that the presence of the flat plate radically 
altered the local heat transfer coefficient around the cylinder especially 
in the classical subcritical Reynolds number range (Reynolds number less 
than 120,000). In the critical Reynolds number range, the overall trend 
in the curves was similar, but the comparative values of the Nusselt 
number at locations of interest were widely differing. In the region 
from 0 to 120 degrees, the values found in this investigation for Reynolds 
numbers in the range between 53,000 and 153,000 were on the order of those 
found by Meyer for a **free" cylinder at Reynolds numbers above 300,000. 
Field's point of maximum heat transfer occurred at the 180 degree mark 
whereas, this investigation showed the maximum value to occur at the 
point of reattaclimcnt of the turbulent boundary layer. Meyer also found 
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icure 10. Comparison of heat transfer results with those of Field at a Reynolds number of 53,000 
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Fiqure 11. Con^^arlson of heat transfer results obtained at a Reynolds number of 153,000 
v/.it]i those of other investiaators . 
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Figure 12. Average Musselt number versus Reynolds number. 
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the maximum value to occur at the reattachment point, but this phenomenon 
was also observed at Reynolds numbers aoove'"'!St;U/000 . 

It is also interesting to note that, despite the lack of correlation 
between this data and Field's data, the plots of average Nusselt number 
versus Reynolds number have nearly the same slope in the Reynolds number 
range between 50,000 and 150,000, with Field's average values being 
higher than those obtained in this investigation at all Reynolds numbers. 

B. PRESSURE EXPERIMENT 

Raw data and sample calculations from the two pressure runs are con- 
tained in Appendix B. This section contains graphical results of the 
pressure coefficient versus angular location for Reynolds numbers of 
105,000 and 156,000 (Figure 13) , and a comparison of the pressure 
coefficient curve at a Reynolds number of 156,000 with rhose obtained 
by Meyer [3] at Reynolds numbers of 153,000 and 495,000 (Figure 14). 

The pressure distributions at the two Reynolds numbers tested were in 
good agreement with the heat transfer results obtained at similar Reynolds 
numbers. This is especially evident in the higher Reynolds mamber case. 

For this case, the pressure coefficient reached a minimum at 73 degrees, 
and the subsequent upward slope contained a double inflection character- 
istic of the separation "bubble” region between 75 and 85 degrees. The 
first inflection point occurred at approximately 78 degrees v;hich corre- 
sponds precisely to the laminar separation point found in the heat transfer 
experiment. Likewise, the second inflection point occurred at approximately 
90 degrees which corresponds closely with the 92-dogree location of the 
reattachment point found in the heat transfer experiment. Additionally, 
the pressure coefficient curve b(?gan to flatten at approximately 135 
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Figure 14. Con.parison of Pressure results at a Reynolds number of 156,000 
with those of Meyer at Reynolds Numbers of 153,000 and 495,000 



/ 



degrees which corresponds favora±)ly with the point of final separation 
observed at 140 degrees in the heat ^ Finally, fluctua- 

tions in the pressure coefficient curves were noted in the regions cor- 

f 

responding to the locations of the "trapped" vortices which reinforces 
the heat transfer data in this region. 

A comparison of the pressure coefficient curves obtained at a Reynolds 
number of 156,000 with those of Meyer obtained on one side of a "free" 
cylinder at Reynolds numbers of 153,000 and 495,000 indicated that the 
behavior exhibited by the cylinder attached to the flat plate more closely 
patterned that of a "free" cylinder at much higher Reynolds numbers. In 
the region between 0 and 180 degrees, the magnitude of the pressure 
coefficients obtained at a Reynolds number of 156,000 were nearly identical 
to those of Meyer obtained at a Reynolds number of 495,000, but offset 
approximately 10 degrees throughout the entire region. This correlates 
with 'the fact that the magnitudes of the heat transfer coefficient at this 
Reynolds number followed closely those of Meyer at Reynolds numbers above 
, 300,000, 

C. FLOW VISUALIZATION EXPERIMENT 

Because of the limitation in the pump capacity, the highest Reynolus 
number which could be obtained in this experiment was approximately 10,000. 
This was v;ell bolov/ the Reynolds number range examined in both the heat 
transfer and pressure pxperiments, and, as a result, correlation was 
difficult. However, the behavior observed agreed well with the proposed 
theory for the subcritical range. 

Referring to the regions of the cylinder described in the proposed 
theory. Region I behaved aJinost exactJy as predicted. The flov/ was 
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observed to separate from the flat plate and reattach to the cylinder at 
a point just aft of the 0 degree inarL. separated region, a 

large turbulent separation vortex was observed, and a series of small, 
circular, three-dimensional vortices were clearly evident at approximately 
30 to 60 degrees from the point of attachment of the cylinder to the 
plate. 

Region II exhibited behavior characteristic of subcritical flow around 
a "free** cylinder. No separation **bubble*' was noted, probably due to the 
fact that the Reynolds number was too low. However, laminar boundary layer 
separation was clearly observed at approximately 80 degrees. This behavior 
was as predicted by the proposed theory for subcritical Reynolds numbers. 

Region III was highly turbulent and three-dimensional and was very 
difficult to interpret with certainty. It did appear, however, that a 
separated region with its characteristic separation vortex was present. 

No ** trapped** vortices were observed, but the flow appeared to be somewhat 
stagnant in the region extending approximately 30 degrees forward of the 
point of attachment of the cylinder to the plate. 
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VI. CONCLUSIONS 

The theory for the flow around a right circular cylinder attached to 
a flat plate in a crossflow of air which was proposed in this investigation 
appears to accurately describe the actual conditions existing at Reynolds 
numbers based on cylinder diameter from approximately 53,000 to 153,000. 

The heat transfer data obtained correlates well with the pressure data 
obtained, and both are in agreement with the proposed theory in a qualita- 
tive sense. Quantitatively, both the heat transfer data and the pressure 
data for Reynolds numbers in the range of 53,000 to 153,000 correlates in 
the region between 0 and 180 degrees with the data obtained by Meyer [3] 
for flow around a "free** right circular cylinder in a crossflow of air at 
Reynolds numbers in the range 300,000 to 500,000. 

The flow visualization exhibited consistency with the proposed theory 
for subcritical flow, but was inconsistent in some regions with the experi- 
mental results. This was probably due to the fact that the flow visualiza- 
tion experiment was conducted at a much lower Reynolds number than were the 
heat transfer and pressure experiments. It is recommended that an attempt 
be made to refine the existing flow visualization apparatus such that higher 
Reynolds number flows can be attained. This would give a clearer picture 
of the flow at Reynolds numbers in the range of 50,000 to 150,000, and 
correlation with the other data could be made more meaningful. 

The liquid crystal technique developed by Field (2] again proved to 
be an excellent method for quickly and easily obtaining heat transfer data 
in forced convection environments. Since the problem addressed in this 
work had not been investigated previously, no measure of comparison is 
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available for the values of the local Nusselt number obtained. But^ the 
qualitative information displa^JS^^ r:.;r:; 4 Stals was verified 

both by the pressure data and the flow visualizatxon. Continued experiments 
involving liquid crystals are highly encouraged. 

The results obtained in this investigation suggest several areas in 
which follow-up type investigations might be conducted. First, the present 
investigation should be extended over a greater Reynolds number range with 
primary attention being focused on Reynolds numbers below 50,000. This is 
considered necessary in order to ascertain at what Reynolds number the flow 
transitions from siobcritical to critical for this particular flow configura- 
tion. It is felt that such a point does exist because the results of the 
flow visualization conducted in this investigation suggest sxibcritical 
behavior at a Reynolds number of 10,000. 

Secondly, an effort should be made to determine at what distance away 
from the flat plate the cylinder can be moved and not feel the effects of 
the plate, i.e. at what distance away from the plate will the cylinder 
behave as a "free" cylinder. 

Finally, a study could be undertaken with the same configuration, but 
with attention being focused on the flat plate rather than the cylinder. 

' This study could be made in the same light as those conducted on forward- 
facing and rear- facing steps. 
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. A . . . . 

HEAT TRANSFER DATA SNU 'UAra ‘RtOoCnON 

The values of the local Nusselt number at various angular locations 
were determined using the relationship between the convective heat trans- 
fer coefficient, h , and the surface heat flux produced by Joulean heating, 

C ' 

V^/RA . This relationship was developed by Meyer [3] by performing a 
s 

simple energy balance on an elemental control volume of the Temsheet. 

The result is 

vVra + (kt/r ^)d^T/d0^ 
s o 

h = h 

^ fT-T„) 



where 

h = convective heat transfer coefficient 
c 

V = voltage applied to the test section 
R = electrical resistance of the test section 
A = surface area of the test section 

k = thermal conductivity of the test section material 
t = thickness of the test section 
r^ = radius of the cylinder 
T = temperature at angular location 0 
T = temperature of the air freestream 

oo 

0 = circumferential location on the cylinder 

2 2 

h = radiation heat transfer coefficient = aF (T+T ) (T +T ) 
r l-- z 

v;here 

0 - Stef an-Boltzman constant = .1714 x 10 BTU/hr-ft - R 

2 “ radiation exchange factor bctv;een cylinder and surroundings 
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Once h has been determined/ the Nusselt number can be obtained from 
c 

the equation ‘ 



Nu = 2h r /k . 

c o air 



The Reynolds number is 

Re = U D/v 

CO 

where 

U = velocity of the air freestream 
00 

D = diameter of the cylinder 

V = kinematic viscosity of air 

Field [2] determined that the product kt for Temsheet was on the order 
of 0.002 BTU/hr - °F; hence, it was considered negligibly small for purposes 
of this investigation. 

A sample calculation of the Nusselt number is provided below using the 
data obtained at a Reynolds nmber of 153,520 and an angular location of 
330 degrees. 

^ 1-2 = 

a = .1714 X 10~® BTU/hr-ft^-°R'^ 

T = 111.0°F = 571°R 

T = 73.1°F = 533. 1°R 

00 

U =76 ft/sec 

GO 

V = 37.4 volts 

R =13.1 ohms 

D/v = 2020 sec/ft 

2r /k . =22.0 hr-°F-f t^/BTU 

o air 

A = 0.528 ft^ 
s 

Re = U D/v = (76 ft/sec) (2020 soc/ft) = 153,520 

u> 

h = h + h = ^ 

^ ^ RA (T-T ) (13.1) (.528 f ) ( 37 . 9°F) ( . 293 watts/BTU/hr ) 

s 
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h =18.2 BTU/hr-ft^-°F 

= (.9) (.1714 X 10”®BTUAr-ft^-'^R’y(briV5i5i.l) (571^+533.1^) °R^ 

h =1.04 BTU/hr-ft^-°F 
r 

h =17.17 BTU/hr-ft^-°F 
c 

Nu = (17.17BTU/hr-ft^-°F) (22.0 hr-ft^-°F/BTU) 

Nu = 377.7 

The data sheets for all nine heat transfer runs follow. 
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DATA 
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DATA 
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DATA 
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APPENDIX h 



PRESSURE DATA AND DATA REDUCTION 

The values of the pressure coefficient were determined from the 
formula 

P- - P 
0 s 



^ pu 



where 

C = pressure coefficient 
P 

= static pressure on the cylinder surface at angular location 0 

0 

= static pressure in the wind tunnel test section 
p = density of air 

= air free stream velocity 

2 

Meyer [3] had determined previously that the quantity ^ pU^ could 
be represented by the difference in static pressure between the wind 
tunnel plenum chamber and the wind tunnel test section. Therefore, 

^ PU ^ = P - p 

OO pc s 

v/here 

P^^^ = static pressure in the wind tunnel plenum chamber 

A sample calculation of the pressure coefficient is provided belov; 

using data obtained at a Reynolds number of 156,180 and an angular location 

of 90 degrees 

P^ = 12.05 in- H^O 

P “4.02 in. H„0 
pc 2 

P = 6.58 in. H^O 
s 2 
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C - 6.58) 

P 



C = -2.14 
P 



The data for the two pressure runs follow. 
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PRESSURE DATA 



Air Speeds j,SO CIH.H2O 
= 51 ft/sec 



Air Temp. = 6i,5 
D//= 30C5 sec/ft 
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PRESSURE DATA 



Air Speeds 3^^ cm.H20 
= 74 ft/sec 



Air Temp.s 62 ®P 
D/V= A-So sec/ft 
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M*PENuia l.' 



UNCERTAINTY ANALYSIS 

Uncertainty calculations were performed to ascertain the degree of 
uncertainty in the final results. The method of Kline and McClintock 

described in reference 19 was used. 

Uncertainties existed in the following measured variables: voltage, 
resistance, cylinder geometry, surface temperature, Temsheet emissivity, 
angular location, manometer pressure drop, all pressures in the pressure 
experiment, and properties of air. The uncertainties associated with the 
latter were considered negligible. 



TOTAL HEAT TRANSFER COEFFICIENT; h 

The total heat transfer coefficient was calculated from the equation 



h = 



2 

V /R 



A (T-T^) 
s ~ 



Therefore, the uncertainty in the heat transfer coefficient is 
expressed as 



k 



/^r ^ ^ 



/O) 









AT 



For the experimental run at a Reynolds number of 153,520 and an angular 
location of 330 degrees, 

V = 37.4 ±0.1 volts (20 to 1) 



R = 13 . 1 ± 0 . 2 ohms 
AT = 37.9 ± 0.5 °F 
A = .528 ± 0.010 ft" 



(20 to 1) 
(20 to 1) 
(20 to 1) 
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• 027 




h 

h = 18.2 ± ,5 BTU/hr-ft"-^F (20 to 1) 

This result does not take into account the uncertainty in angular loca-- 
tion in this region. This uncertainty is estimated to be ± 3 degrees 
(20 to 1). 

The uncertainty in the total heat transfer coefficient remained rela- 
tively constant around the entire circumference of the cylinder. However, 
the uncertainty in angular location varied from a maximum of ± 5 degrees 
(20 to 1) in the stagnation region and the wake region to a minimum of- ± 1 
degree (20 to 1) in the regions near the laminar separation point and the 
reattachment point. 



RADIATIVE HEAT TRANSFER COEFFICIENT, h 
^ 

The radiative heat transfer coefficient was determined from the 
equation 



2 2 

h = aF, ^ (T + T ) (T + T ) 
r 1-2 ® ® 

Assuming the uncertainty in the temperature measurement is the same 
for both T and T^, and defining = % (T + T^) , 

h = 4 oF^ _T ^ 
r 1-2 m 

Since o is a universal constant, it does not enter into the uncertaint]. 
Therefore, 




The radiation exchange factor, F^ can be shov;n to be equal to the 
emissivity of the Temsheet. Since this value v/cis )iot known, a value of 
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0*9 was chosen based upon the emissivities of comparable materials. The 
uncertainty in this value was assumed to be ± .05 (20 to 1). Hence, at a 
Reynolds number of 153,520 and an 'v of 330 degrees 

-r^ = 0.06 

h 

r 

\ = 1.04 ± 0.06 BTU/hr-ft^-°F (20 to 1) 



NUSSELT NUMBER, Nu 

The Nusselt number was calculated from the equation 

Nu = h D/k = (h-h )D/k 
c r 

Therefore 






Nu 



Nu 



• ft) ■ ft • ft 



Since the uncertainties in the properties of air were considered 



neglible and the uncertainty in the 
small. 




Nu =378 ± 11 



diameter measurement was negliglibly 



2 2 




(20 to 1) 



REYNOLDS NUI-SER, Re 

The Reynolds number was calculated from the equation 

Re = U D/v 

CO 
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Therefore 






The quantity is directly proportional to the square root of the 
wind tunnel micromanometer reading AH; hence 



(X) 0). 

U AH 

00 

= 0.5 



U 



AH 



Again, since the uncertainty in the kinematic viscosity and in the 
cylinder are negligible. 



(l) (x) . 

“ . 0.5 ^ 



Re 



AH 



At a Reynolds number of 153,520 (AH = 3.58 cm. the uncertainty 

in the micromanometer reading was estimated to be ± 0.03 cm. Hence, 



0) 



Re 



Re 



= 0.5(0.03/3.58) 



(X) 



Re 



Re 



= 0.004 



Therefore, 



Re = 153,520 ± 614 



(20 to 1) 



PRESSURE COEFFICIENT, C 

; P 

The pressure coefficient was determined from the equation 



P. - P 

c ^ 

p p - p 
pc S 



AP, 



AP 



pc 
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Therefore 




Since both the pressure differences were measured on the same man- 
ometer bank/ they shared the same uncertainty which was considered to be 
± 0.05 in. H^O (20 to 1). Therefore, at a Reynolds number of 156,180 and 
an angular location of 90 degrees , 




0.022 



C = -2.14 ± 0.05 
P 



(20 to 1) 
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